Orbital Resonance Mode in Superconducting Iron Pnictides 
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We show that the fluctuations associated with ferro orbital order in the dxz and dyz orbitals can 
develop a sharp resonance mode in the superconducting state with a nodeless gap on the Fermi 
surface. This orbital resonance mode appears below the particle-hole continuum and is analogous 
to the magnetic resonance mode found in various unconventional superconductors. If the pairing 
symmetry is s±, a dynamical coupling between the orbital ordering and the d-wave subdominant 
pairing channels is present by symmetry. Therefore the nature of the resonance mode depends on the 
relative strengths of the fluctuations in these two channels, which could vary significantly for different 
families of the iron based superconductors. The application of our theory to a recent observation of 
a new (5-function-like peak in the Big Raman spectrum of Bao.6Ko.4Fe2As2 is discussed. 
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Introduction - For high-temperature superconductors, 
cuprates and iron pnictides, resolving the nature of the 
fluctuations in both the normal and superconducting 
states remains a crucial question as it holds the key to the 
pairing mechanism. In cuprates, due to the strong an- 
tiferromagnetism in the parent compounds, it is widely- 
believed that the antiferromagnetic spin fluctuations are 
the most important ingredient in the pairing mechanism. 
One of the marquee indicators of this is the magnetic res- 
onance mode[H observed in the superconducting state of 
every cuprate. From the BCS theory, the spin-flip sus- 
ceptibility of an electron scattered from fc to A: -I- g in 
the superconducting state gains an extra coherence fac- 
tor which is proportional to (1 — sgn(A(A;))sgn(A(fc-|-5))). 
Since the gap symmetry of the cuprates is c?-wave signi- 
fied by sgn(A(fc)) = — sgn(A(A: -I- Q)), the spin excita- 
tions near q~Q are compatible with superconductivity. 
It can be further shown [2-l6| that a sharp (5- function-like 
resonance mode in the spin-flip susceptibility requires an 
antiferromagnetic spin interaction to pull the resonance 
mode to an energy below the particle-hole continuum. 
In other words, despite some dependence on the detailed 
electronic structure, the existence of the magnetic res- 
onance mode is predominantly determined by the gap 
symmetry and the nature of the spin interaction. As a re- 
sult, it has been identified as an unambiguous [l[ indicator 
that antiferromagnetic spin fluctuations remain strong in 
the superconducting state and consequently might be the 
driving mechanism for superconductivity. 

In iron pnictides, the structural phase transition fol- 
lowed by stripe-like antiferromagnetism is a robust fea- 
ture in the phase diagram. The onset of the structural 
phase transition breaks the C4 symmetry down to C2 
which enables magnetic order at wavevector (tt, 0) and 
orbital order at zero wavevector characterized by un- 
equal occupation in d^z and dyz orbitals. These two 
unique consequences suggest that the dominant fluctu- 
ations might be either orbital-based [7 13] or spin-based 
Indeed, a magnetic resonance mode at wavevec- 
tor (tt, 0) has been observed in the pnictides and has 



been attributed to spin fluctuations Uj]. However, re- 
cent Raman scattering measurements [l9|, a zero wave- 
vector probe, has observed a superconductivity-induced 
peak in the Big channel. Since this peak occurs at zero 
wavevector, it cannot be attributed to spin fluctuations. 
It can, however, arise from a sub-dominant d-wave pair- 
ing channel lol- 21 1. Our key point here is that d-wave 
sub-dominant pairing and ferro-orbital order both have 
the same space group symmetry of the Raman Big mode. 
Consequently, the observed Big mode in Raman could 
arise from orbital fluctuations. We term this mode an 
orbital resonance mode. We work out the details of this 
scenario and show that the observed Big mode could of- 
fer an unprecedented fingerprint of ferro-orbital order in 
the pnictides. The Raman experiments then supplement 
the current indirect evidence in the non -sup erconducting 
states that point-contact spectroscopy 221424 1 and neu- 



tron scattering measurements 2^ 2^ provide for orbital 
fluctuations. 

Criteria for the orbital resonance mode - First we give 
a general discussion for the existence of the orbital res- 
onance mode as was done for the magnetic resonance 
mode. Given that orbital fluctuations are in the charge 
channel, the extra coherence factor from BCS theory now 
becomes ~ (1 -|- sgn(A(fc))sgn(A(fc -I- q))). We imme- 
diately see that the sharp difference from the magnetic 
counterpart is that the prerequisite condition for the or- 
bital resonance mode is sgn(A(A:)) = sgn{A{k+(f)). Since 
we focus on the case with wavevectors near zero {q ^ 0), 
such a condition generally holds in any gap symmetry. In 
the following, we show that the orbital resonance mode 
generally exists provided an effective attractive interac- 
tion is present in the orbital ordering channel. 

Before moving to the microscopic calculations, we 
prove that the orbital ordering and d-wave sub-dominant 
pairing channels are coupled in an s-wave superconduct- 
ing state. This state of affairs obtains for two crucial rea- 
sons. First, because the singlet Cooper pair is a mixture 
of electrons and holes with different spins, the particle- 
hole excitations arc intrinsically coupled to the particle- 
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particle excitations as long as the space group symmetry 
allows. Second, both orbital ordering and the d-wave 
pairing channels have a sign change under a rotation of 
7r/2 but the s-wave ground state does not. Therefore a 
Berry-phase coupling between them is allowed by sym- 
metry as discussed in Ref. This is the key physical 
principle underlying our work here, and we will show that 
such a Berry phase coupling enriches the physics of the 
Big resonance mode. 

Formalism - We start from a general two-orbital model 
of the superconducting state for iron pnictides. The 
model Hamiltonian is given by 



H 
Ho 



Hq + Hsc + Hi, 

'^i'l^ e+{k)i + e-{k)T3 + ea;y{k)Tx 
V^t Dr (t>t , 



Hsc 



k,(T 

E 

fe 



A{k) 



(1) 



where ibt = (c?t dl ) are the electron creation 

k.<T k.xz.a k^yz^a 

operators for orbital xz or yz with spin a, Dr, = 



diag.{E'2 ,Ei ) the eigenenergy, and </)t = (at ) 
are the corresponding eigenvectors. We adopt the model 
proposed by Raghu et. aL[27] for Hq and use the same 
tight-binding parameters. Hsc describes the mean-field 
pairing potential in the diagonalized basis and we as- 
sume the pairing symmetry to be s±, which in the two- 
orbital model is A{k) = Aq cos(fca:) • cos(fcy). Without 
loss of generality, we consider two types of interactions 
Hi = Hoc + Hdsc, 

^OO = ^J2^Ooiq)0{-q)diq), 



At - 



Hdsc = 



(2) 



where 



0{g) 

A\q) 



k.CT 



dk 



cos k^ — cos k 



y 



(3) 



Orbital order is then characterized by {0{q = 0)) ^ 0, 
and Hqo is the effective interaction in the orbital or- 
dering channel. The orbital ordering instability requires 
VooiQ = 0) to be negative and smaller than a critical 
value. Finally, Hdsc represents the interaction for the 
sub-dominant pairing interaction in the d-wave channel. 
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FIG. 1: The dynamical couplings between the orbital ordering 
and d-wave pairing channels. The magnitude of the gap is 
Ao — 0.01, and the edge of the particle-hole continuum 
is 0.0127. A typical behavior of the Berry phase coupling: 
Xi^) ~ ^ can be seen. 



which is generally believed 2^31| to exist in iron-based 
superconductors due to their unique Fermi surface topol- 
ogy. Because we take the ground state to be the s± su- 
perconducting state with which the d-wave pairing chan- 
nel is competing, the instability toward d-wave supercon- 
ductivity will occur only if Vd^q = 0) exceeds a critical 
value. In other words, the ground state of s± supercon- 
ductivity is stable as long as Voo{q = 0) > —Vqq and 
Vdiq = 0) > -V^, where VSo,V^^ > 0. Note that both 
terms can be obtained from a microscopic multi-orbital 
Hubbard Hamiltonian. 

Since our main purpose is to study the behavior of 
the fluctuations due to Hqo and Hdsc in a stable s± 
superconductor, we choose Voo{<i — 0) and Vd[q = 0) 
to be close but not exceeding critical values. The pro- 
cedure for analyzing our model Hamiltonian is well es- 
tablished. First, the Bogoliubov transformation is per- 
formed to diagonalize {Hq -\- Hsc)- Second, the random- 
phase approximation is applied to the two-particle corre- 
lation functions to express them in terms of the Bogoli- 
ubov quasiparticles. The relevant susceptibility is given 
by 



(4) 



X°(g, w), w), and are all 3 x 3 matrices. The bare 
correlation functions, x^{q,uj), are defined by 



(5) 



where 



Aiit) 



J{Ho+Hsc)t/h 



0{q)e 



-i{Ho+Hsc)t/h 



(6) 
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A2{t) = e''(-f^''+^s^)*/''(A''(-q) - A''1'(g))e-*(^«+^^c)tA^ 



and Uq- is the effective interaction kernel, 



(7f = 



Vooiq) 

n YiM. 
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(7) 



To assist with the analysis, we divide the d-wave com- 
peting pairing channel into phase {A2) and amplitude 
(A3) modes in order to enforce the time reversal symme- 
try in each of the of Ai channels. It will be shown later 
that because all of these three modes couple to each other 
due to the nature of the superconductivity and the space 
group symmetry, any theory ignoring these couplings is 
incomplete and might overlook important physics. For 
the details, please see the Supplementary Material. 

Results - From now on, we study the case of zero 
wavevector {q = 0) and neglect the subscript q. In princi- 
ple, our formalism can be applied for finite q, but detailed 
knowledge in the (f-dependent interaction will be neces- 
sary in order to obtain the correct dispersion of the res- 
onance mode. For proof of principle, an analysis of the 
zero- wave vector case suffices. It is now clear, from the 
form of [x°(g, w)] that the orbital ordering and the 

d-wave pairing channels are coupled because [x'''^)] ^ is 
non-zero even for i ^ j. In Fig. [TJ we plot [x*'(a;)] and 
[x''(w)] for u! smaller than the particle-hole continuum 
edge « 2A{kp), which show the typical behavior for 
a Berry phase coupling: x{^) 6]. A similar coupling 
occurs in the magnetic resonance mode of the cuprates 
in which the 7r-particle and spin-fiip channels are coupled 
by symmetryUSQ- 

Since an s± superconductor is fully gapped on the 
Fermi surface, any of the two-particle correlation func- 
tions should have zero spectral weight below flph unless 
there exists a resonance mode with frequency iOres which 
is a solution to 



Det([xKe.)] ') =0. 



If Imx''(w) has a sudden increase from zero at a; = ilph, 
Rex°(w) develops a weak logarithmic divergence at w = 
ilpft, as can be shown from the Kramers-Kronig relations. 
In this case, Eq. ([8]) is satisfied if there is an effective 
attractive interaction in the orbital ordering channel, i.e., 
^00 (0) < 0. In other words, an orbital resonance mode 
can generally arise if the fluctuations from the orbital- 
ordering transition persist as the transition to the fully- 
gapped superconducting state obtains. This orbital res- 
onance mode is an analogue of the magnetic resonance 
mode in the cuprates. 

In Fig. [21 we plot the spectral weight of the bare cor- 
relation function in the orbital-ordering channel. Clearly 
a jump from zero to a finite value occurs at the edge of 
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FIG. 2: The bare correlation function for the orbital-ordering 
channel. The magnitude of the gap is Ao = 0.01, and the edge 
of the particle-hole continuum Qph is indicated by the dashed 
line. To obtain a smoother plot, we included a broadening 
factor S — 0.0002 in our computation of the imaginary part 
of the bare correlation function. 



the particle-hole continuum. This signifies that oj^es Sii- 
ways exists below flph for Vbo(O) < 0. Let us discuss 
first the case without a residual c?-wave pairing interac- 
tion Vd(0) = 0. Generally speaking, if the interaction 
strength |Vbo(0)| is small, ojres/^ph ~ 1 and the spec- 
tral weight in this resonance mode is small. On the other 
hand, there exists a critical strength of |Vbo{0)| for which 
Eq. [8] yields a solution of w^es — (in this calculation, 
this critical strength is |Vbo(0)| ~ 1-2). This marks the 
instability towards the orbital ordering transition. Since 
we are only interested in the superconducting state with- 
out orbital ordering, we have < ujres/^ph < 1- This 
behavior is confirmed by our RPA calculations which are 
summarized in Fig. [3] 

Note that the orbital resonance mode discussed so far 
will not be present for a superconducting gap with nodes. 
The reason is that the orbital resonance mode emerges 
near zero wavevector. Consequently, the particle-hole 
(^8) continuum edge is roughly 2Amm, where Amin is the 
minimal gap on the Fermi surface. If Amin — 0, which 
is the case of the gap symmetry with nodes on the Fermi 
surface, Eq. [8] can never be satisfied. Instead, a peak 
could still arise, but it will be damped or even com- 
pletely washed out by the particle-hole continuum, the 
precise details of which depend on the electronic struc- 
ture of the material and thus is not symmetry protected. 
The aforementioned case might be realized in cuprates 
in which nematic fluctuations (identical to the orbital 
fluctuations [2j|) are argued to be present [32, |33| and the 
d-wave gap symmetry has nodes on the Fermi surface. 
This is in a sharp contrast to the magnetic resonance 
mode which usually appears near flnite wavevector, for 
example, (tt, tt) for the cuprates and (tt, 0) for iron pnic- 
tides. In this case, the particle- hole edge, flph, could be 
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finite for botii nodal and nodeless gap symmetry, and 
the existence of the magnetic resonance modes is deter- 
mined primarily by the sign difference between A{k) and 
A{k + q). 

The Berry phase coupling gives rise to rich physics in 
the context of the resonance mode. If only one of these 
two channels has strong fluctuations, the resonance mode 
has an increased weight in the dominant channel. There- 
fore, the Big mode can be viewed as either an orbital 
resonance mode or a c?-wave excitonic pairing mode de- 
spite the non-zero dynamical coupling between them. If 
both channels, however, are equally strong, then the res- 
onance mode has comparable weight in both channels. 
As a result, the nature of the resonance mode becomes 
plasmonic[l,0|, and the resonance frequency tends to be 
smaller as shown in Fig. [3l It is interesting to point out 
that in iron pnictides, the relative strength of the orbital 
ordering and the d-wave pairing fluctuations could vary 
significantly from material to material. This means that 
the nature of the resonance mode discussed above could 
be different for different families of the iron based super- 
conductors. In light of this observation, we argue that 
the interpretation of a (5-like peak recently found in the 
Raman scattering measurement on Bao.6Ko.4Fe2As2 and 
Rbo.8Fe1.6Se2 needs further consideration. 

While Kretzschmar et. al. [l9j] have given a detailed fit 
to their data based on a model consisting of only com- 
peting pairing potentials, we would like to offer an alter- 
native viewpoint based on the present theory. We note 
that from Fig. 3 in Ref.[l^, all of the peaks observed in 
Rbo.8Fc1.6Se2 appear both above and below the super- 
conducting transition temperature Tc, while in Fig. 4(c), 
the peaks in the Big symmetry in Bao.6Ko.4Fe2As2 are 
present only in the superconducting state. This indicates 
that only Bao.6Ko.4Fe2As2 has truly superconductivity- 
induced collective modes of Big symmetry. Since the or- 
bital order operator, 0{q = 0), defined in Eq. [3] changes 
sign upon a rotation by n/2, it has Big symmetry and 
consequently the Big Raman spectroscopy is the ideal 
probe of the orbital resonance mode discussed in this pa- 
per. The first peak at w = 190cm^^ is the feature of 
the particle-hole continuum edge, and we interpret the 
second peak at lu — 140cm~^ to be the orbital reso- 
nance mode. For Rbo.8Fe1.6Se2, because a structural 
phase transition is not observed and theoretical calcu- 
lations also suggest an absence of ferro-orbital ordering 
3J-|36[, the orbital resonance mode should not exist in 
this material. 

Conclusion - We have shown that a resonance mode 
arising from fluctuations of the ferro-orbital ordering 
channel can exist in the superconducting state with s± 
gap symmetry. This mode has a dynamical coupling to 
the d-wave sub-dominant pairing channel due to sym- 
metry. Consequently, the nature of the resonance mode 
depends crucially on the relative strength of the fiuctua- 




FIG. 3: The ratio of the resonance frequency to the particle- 
hole continuum edge uJres/Qph as a function of |Voo(0)| 
for different values of Vd(0). At small interaction strength, 
ujres/^ph ~ 1. The instability of orbital ordering in the su- 
perconducting state occurs at Voo(O) « —1.2, which has the 
solution of Ures = (not shown). Larger |Vd(0)j tends to 
shift the resonance frequency to smaller values. Since we are 
only interested in the superconducting state without orbital 
ordering, < ujres/^ph < 1- 



tions in these two channels, which could vary significantly 
among the members of the iron-based superconductor 
family. If our scenario is correct, then Raman experi- 
ments on pnictides which have no structural transition 
should also not display the Big resonance mode. Hence, 
our theory is inherently falsifiable. 
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Supplementary Materials 

Here we provide all the details of the formalism used 
in the paper. First, the tight-binding Hamiltonian Hq is 
defined as 



= y^di bt (j>j: 



k,(T 



(9) 



k.a 



= (di ,di ), br = diag.(Ee ,e1 ), and 
(bL = (at , Bl ). We introduce a unitary transforma- 
tion JJr such that 



Ul^-D., (10) 



and it can be derived that 



where 



^ f cos 9k sin Ok 
^'<^ \ - sin^fe cosOk 



cos29k = — , sm26'fc = — 

H{k) H{k) 



(11) 



(12) 



The eigenvalues are 

^ e+{k)+H{k) , El^e+{k)^H{k), 



{k)+ey,{k) 



e+{k) EE 



H{k) EE ^e^{ky+ely, (13) 
and the corresponding eigenvectors can be expressed as 

The next step is to perform a Bogoliubov transforma- 
tion for H' = Hq + Hsc- We then define 



I3,k 

Then we can rewrite 
Ho +Hsc = y] 



E^f3 + A^{k)fi 



5' 



where 



k,^—a,fi 
k,fi—a,l3 



a.k - ^ +.k ' 

$t , EE (/3^^-.t,/3«^_t) 

I3,k ^'^+,k ' ^-,k ' 



(16) 



(17) 
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are the Bogoliubov quasiparticlcs in the a or /3 bands. 
The energies and the corresponding eigenvectors of the 
BogoUubov quasiparticles are 



where 



k 



sin 2uj t 



A^(fc) 



(20) 



fe 



(18) 



where wc have introduced another unitary transforma- 



tion U^'^ such that 



Ilk 



COSCJ t 



(19) 



All that is left is a series of long but straightforward 
caluclations. Wc express all the three channels Ai,A2, A3 
in terms of the Bogoliubov quasiparticles {a^,(3^'~'}, 
and then the bare correlation functions are 



AtA 



lAJ 



+ 



Bi Bj 



Ct^^ tvjj + i5 + {E'^ci^ - q) + E'^ci'^)) huj + iS-{E^cik-<l) + E'^cik))- 

k 



+ [- 



CiCj 



+ 



DiDj 



nw + i5+ {E^scik -q) + Elcik)) fkj + i5- (^f ^(fc - q) + E^^ik)) 



-'sc 

1 



+ {EiEjx\ 5 \ 3 1 

huj + i5+{E'^c{k-q)+E^sc{k)) + id - {E^cik - q) + E^ciW 

+ FiFjX \ ^ 1 1 (21) 

'hw + i6-{E^scik-Q} + E^cik)) huj + id + (Elcik - q) + E^cik)) 



where 
Ai -- 
A2 -- 
A3 = 
Ci = 
C2 -- 



Bi = cos{ek-q + Ok) X sin (w^ j + 

-B2 = -d^ X cos (w^ g + u}^^k-q) 

— B3 = — dr: X cos ( W r — W r -1 

Di = - cos(6ife_, + 9k) X sin (u^ j^ + 

-D2 = -dg X cos (u^ j; + 



C3 = 

El = 

E2 = 

Fi = 
F2 



--D3 = -d^ X cos (w^ g - 
sin(6lfc_g + Ok) X sin (^w^ g + 

^3 = 

sin(6ifc_, + Ok) X sin {uj^j^ + uj^^k-g) 



= F3 = 0. 



(22) 



